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Operation of a Superconducting Quantum Interference Device (SQUID) made of stacked Josephson junctions 
is analyzed numerically for a variety of junction parameters. Due to a magnetic coupling of junctions in the stack, 
such a SQUID has certain advantages as compared to an uncoupled multi-junction SQUID. Namely, metastability 
of current-flux modulation can be reduced and a voltage-flux modulation can be improved if junctions in the 
stack are phase-locked. Optimum operation of the SQUID is expected for moderately long, strongly coupled 
stacked Josephson junctions. A possibility of making a stacked Josephson junction SQUID based on intrinsic 
Josephson junctions in high-T c superconductor is discussed. 
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1. Introduction 

Layered high-T c superconductors (HTSC) rep- 
resent natural stacks of atomic scale intrinsic 
Josephson junctions (IJJ's) 0. Indeed, an in- 
tcrlayer spacing of 15.5 A was estimated from a 
periodic Fraunhofer modulation of Fiske steps in 
Bi-2212 mesas [||. At present, high quality SIS- 
type IJJ's can be fabricated using HTSC mesa 
structures J3|||]. IJJ's with their record large I c R n 
values ~ lOmU are attractive candidates for cry- 
oelectronics. 

Here I analyze a possibility of making a SQUID, 
based on IJJ's. Because the interlayer spacing in 
HTSC is very small, it is likely that such a SQUID 
will contain several stacked Josephson junctions 
(SJJ's), as shown schematically in Fig. 1. From 
previous studies it is known that operation of a 
multi-junction SQUID, consisting of series con- 
nected uncoupled junctions, is problematic due 
to the existence of metastable states and a small 
current-flux modulation The main differ- 

ence between series connected and stacked junc- 
tions is that SJJ's are coupled with each other, 
which may result in qualitatively different behav- 
ior of SJJ's, see eg. recent review g. 

The main question which is raised in this paper, 
is whether the operation of a multi-junction 



SQUID can be improved by coupling and phase- 
locking in stacked Josephson junctions. 

2. Magnetic coupling in stacked junctions 

I will consider a magnetically (inductively) cou- 
pled stack of N junctions. The coupling is 
achieved via common electrodes and is strong 
when they are thinner than the London penetra- 
tion depth, As (which is certainly the case for 
HTSC IJJ's). Properties of magnetically cou- 
pled SJJ's are described by a coupled sine-Gordon 
equation M, 



cpl = A ■ Jisin((fii) 



(1) 



Here (pi are gauge invariant phase differences, 
'primes' denote in-plane spatial derivatives, A 
is a symmetric tridiagonal N x N matrix with 
nonzero elements: A^j-i = —Si/ 'Aj; A^j = 
Aj/A ; ; = —Si+i/Ai, where A* = 

U + XsiCotti(j£) + Xsi+icoth^j^j, $ = 

Xsicosech (^j^-^j, di and ti are thicknesses of su- 
perconducting electrodes and insulating barriers, 
respectively, and 3b is a bias term [jlO| . 

The coupling strength S is determined by off- 
diagonal elements of A. For identical junctions, 
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S ~2ch- 1 (d/\ s ). 



(2) 
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Figure 1. A schematic view and a diagram of a 
stacked Josephson junction SQUID. 



3. Short junctions, Ls < Xj 

According to Eq.(l), SJJ's are coupled via the 
second derivative of phase. The characteristic 
length for the phase variation is given by the 

Josephson penetration depth, Xj 



8ir 2 J e A s O- 

For short SJJ's, Ls < Xj, magnetic coupling is 
inefficient and a SJJ SQUID is described by sim- 
ilar expressions as an uncoupled multi-junction 
SQUID, see the diagram in Fig. 1. 



i = h + i2\ 

Jiisin(tpu) = ii, J2jsin(if2j) = i - iv, 



J 



(pii — arcsin \-j^sin(tpik) j + '2ttii. 

. ( 'i - Jiksin(tpikY . 

ifi2j = arcsin | + Zirm; 

V J 2j 
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Here subscripts 1 and 2 correspond to left and 
right legs of the SQUID, respectively, <& = <I> e — 
i\L\ + 12^2 is the total magnetic flux in the 
SQUID loop, <E> e is an external flux and Lx 2 are 
inductances of the SQUID legs. All variables are 
expressed in terms of the phase difference in the 
weakest junction " lk" . 

The critical current vs. flux characteristic of 
the SQUID with N = 3 uncoupled (L s = 0) junc- 
tions per leg with equal J c 's is shown in Fig. 2 
by thin lines. It is seen that I c (&) consists of a 
number of lobes. Each lobe has a width iV$o (not 
$o as in a conventional SQUID) . This can be eas- 
ily understood from Eq.(6): due to a summation 
in the right-hand side, the overall flux increases 
N— fold as compared to a conventional SQUID 
with a single junction per leg. Nevertheless, the 
overall periodicity of the I c {®) pattern is one flux 
quantum 4>o- This is caused by the appearance 
of equivalent lobes shifted by $o> corresponding 
to 2tt shift of phase in the junctions, see Eq.(6). 

From such a simple analysis we can see two 
main disadvantages of the uncoupled multi- 
junction SQUID: (i) The SQUID possesses sev- 
eral metastable states and I c (&) is not unique, 
(ii) The overall modulation of I c (&) (envelope of 
all lobes) strongly decreases due to a mismatch 
between the width of the individual lobe (iV$o) 
and the spacing between lobes ($o)- 

4. Long junctions, Ls > Xj 

Magnetic coupling is effective in long stacks. 
However, as far as SQUID application is con- 
cerned, the stack should not be too long otherwise 
flux penetration inside the stack becomes com- 
plicated due to the appearance of multiple quasi- 
equilibrium fluxon modes 1 10 1^^] . The optimum 
stack should have a moderate length, Ls ~ 4A,/, 
so that magnetic coupling is effective, but fluxon 
modes in SJJ's are not pronounced. 

Static current-flux characteristics of a SJJ 
SQUID with the geometry shown in Fig.l were 
simulated numerically using a finite difference 
method with successive interactions. For each 
stack, Eq. (1) was solved together with Eqs.(3,6) 
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Figure 2. Simulated I c ($ e ) dependence of a 
SQUID based on TV = 3 strongly coupled SJ Js for 
different junction lengths. It is seen that metasta- 
bility of / c ($ e ) is reduced in moderately long SJJs 
SQUID due to magnetic coupling of junctions. 



Figure 3. Simulated J c ($ e ) dependence of a 
SQUID based on N = 3 SJJs for a) moderately 
coupled and b) strongly coupled SJJs. It is seen 
that strongly coupled SJJs SQUID becomes more 
stable and less susceptible to the spread in junc- 
tion parameters. 



and with boundary conditions that magnetic in- 
duction is equal to external field on the outer 
edges and to the field inside the SQUID loop for 
the inner edges of SJJ's. Simulations were made 
for N = 3 SJJ's in each SQUID leg. The outer 
electrodes were assumed to be thicker, di t 4 > c?2,3, 
in order to model bulk top and bottom bars of the 
SQUID loop, see Fig. 1. 

First we analyze the dependence of current-flux 
modulation on the in-plane length of the stack, 
Ls- Fig. 2 shows simulated i c ($ e ) patterns for 
a SQUID with strongly coupled SJJ's, S ~ 0.5, 
for different Ls- Parameters of SJJ's are: d\ = 
0.2Aj, d 2 , 3 = O.OlAj, d 4 = OAXj, ii_ 4 = O.OlAj, 
Ag = O.lAj, critical current densities of all junc- 
tions are identical. In Fig. 2 a) J c ($ e ) pat- 
terns are shown for a fixed SQUID loop length, 
L = 20Aj, for short SJJ's Ls = 0.5A,/ (open cir- 
cles) and moderately long SJJ's Ls = 4Aj (solid 
squares). For comparison, I c ($ e ) for the uncou- 



pled multi-junction SQUID (Ls = 0) is shown by 
thin lines. It is seem that I c (&e) of the SQUID 
with short SJJ's exhibit similar metastability as 
the uncoupled multi-junction SQUID due to in- 
efficient magnetic coupling in short SJJ's. How- 
ever, magnetic coupling in moderately long SJJ's 
strongly reduces metastability of I c (& e ). An ad- 
ditional large scale modulation of I c ($ e ) in mod- 
erately long SJJ's SQUID is due to flux quanti- 
zation inside SJJ's. 

In order to eliminate differences caused by the 
finite junction length effect, a short SJJ SQUID 
with a smaller SQUID loop Lq — 2\j was also 
simulated. Results are shown by open circles in 
Fig. 2 b) along with the data for a moderately 
long SJJ SQUID from Fig. 2 a). Now finite 
length effects are similar for both short and long 
SJJ SQUID's. However, metastability is much 
more pronounced for a short SJJ SQUID, while 
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it is strongly suppressed in a moderately long SJ J 
SQUID. 

Next we consider the influence of critical cur- 
rent nonuniformity and the dependence on the 
coupling strength. In Fig. 3 a) and b) current- 
flux characteristics of moderately (S = 0.3) and 
strongly (S = 0.5) coupled SJJ SQUID's, respec- 
tively, are shown for identical J c 's (solid symbols) 
and for a 20 % spread in J c 's (open symbols). In 
Fig. 3 a) magnetic coupling is reduced due to a 
larger thickness of middle electrodes, see Eq.(2), 
— O.lAj = As, other parameters are the same 
as in Fig. 2. In the case of nonuniform junctions, 

J C 2 = l.Uei, J C 3 = 1.2J C \. 

From Fig. 3 it is seen that I c {$ e ) of a moder- 
ately coupled SJJ SQUID is metastable and the 
performance is further deteriorated by the spread 
in J c 's. On the contrary, strongly coupled SJJ 
SQUID is more stable and less susceptible to a 
spread in junction parameters. 

The increased stability of strongly coupled SJJ 
SQUIDs can be understood in terms of phase- 
locking of junctions in the stack. Indeed, the 
energy difference between states is increased ap- 
proximately proportional to the number of phase- 
locked junctions. This was confirmed experimen- 
tally by analyzing the "escape rate" in HTSC 
IJJ's jl2|]. Another important advantage is that 
phase-locked junctions switch simultaneously into 
the resistive state ][2| . This creates a possibil- 
ity for compensation of a reduced current-flux 
modulation in the multi-junction SQUID by an 
increased voltage-flux modulation due to phase 
locking of junctions in the stack. 

5. Conclusions 

In conclusion, operation of a stacked Joseph- 
son junction SQUID is studied numerically. I 
demonstrate that the coupling between stacked 
Josephson junctions can improve operation of a 
multi-junction SQUID: 

(i) Metastability of i" c ($) is reduced in SJJs 
SQUID due to magnetic coupling of junctions. 

(ii) Strongly coupled SJJs SQUID becomes 
more stable and less susceptible to a spread in 
junction parameters. 

(hi) Small I c {&) modulation can be compen- 



sated by an increased voltage-flux response if 
junctions in the stack are phase-locked. 

Finally, some remarks about the optimum 
SQUID design: Magnetic coupling is effective in 
long SJJs. However, SJJ's should not be very 
long, otherwise multiple quasi-equilibrium fluxon 
modes in SJJ's would strongly complicate flux 
penetration into the SQUID loop The 
best performance is expected for moderately long 
SJJ (L s ~ 4Aj - 2 - 3^m for Bi-2212). In order 
to achieve reasonable I c (&) modulation the num- 
ber of SJJ's should be kept small (less than ten). 
This is difficult for HTSC SQUID with the geom- 
etry shown in Fig. 1. The solution would be to 
place the SQUID loop in the afc-plane and form 
stacks of intrinsic Josephson junctions by cutting 
trenches, restricting current flow in the afe-plane. 
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